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We study resonant tunneling through quantum-dot systems 
in the presence of strong Coulomb repulsion and coupling to 
the metallic leads. Motivated by recent experiments we con- 
centrate on (i) a single dot with two energy levels and (ii) a 
double dot with one level in each dot. Each level is twofold 
spin-degenerate. Depending on the level spacing these sys- 
tems are physical realizations of different Kondo-type models. 
Using a real-time diagrammatic formulation we evaluate the 
spectral density and the non- linear conductance. The latter 
shows a novel triple-peak resonant structure. 
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Quantum transport of electrons through discrete en- 
ergy levels in quantum dots has been studied both ex- 
perimentally and theoretically In small dots 
the Coulomb repulsion is strong and may yield the dom- 
inant energy scale, which gives rise to Coulomb block- 
ade phenomena (see e.g. Refs. for single dots and 
Refs. 
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for multiple dots). At high temperatures, or 
when the coupling between the dots and metallic leads 
is weak, the electrons tunnel through the system sequen- 
tially. At low temperatures and for strong coupling to 
the leads resonant tunneling processes contribute signif- 
icantly to the current ||,|l^-|l^ . 

The purpose of this Letter is to investigate electronic 
transport phenomena through two ultrasmall quantum- 
dot systems: (i) a single dot with two energy levels and 
(ii) a double dot with a single level in each dot (see 
Fig. 1^) . We will show that the two models can be mapped 
onto each other. Thus, it is sufficient to concentrate on 
one, the single-dot model with two levels. In both cases 
we assume strong Coulomb repulsion between the elec- 
trons, such that we have to consider only two charge 
states, e.g. the empty or the singly-occupied dot sys- 
tem. We further allow for strong coupling between the 
quantum dot(s) and the metallic leads adjacent to them. 
This gives rise to resonant tunneling and ~ as has been 
shown for a single dot with one energy level |p^|l7| - 
nonequilibrium Kondo effects. In the following we con- 
sider the case with vanishing magnetic field. Accordingly 
the single-electron states are spin-degenerate, which is a 
prerequisite for Kondo effects. We suggest that the mod- 
els could be realized experimentally in a setup resembling 
those in Refs. |l|-|2^. 

A systematic description of resonant tunneling phe- 



nomena has been developed in Refs. ||3|,|15|,^ in terms of 
a real-time diagrammatic technique. It has been applied 
to the problems of electron transport through a single- 
level quantum dot as well as through a small metal island, 
both connected to reservoirs via tunnel junctions. In the 
quantum dot Kondo-type effects manifest themselves as a 
zero-bias maximum in the differential conductance. Fur- 
thermore, a zero-bias minimum has been predicted for 
the case where the level lies above the electrochemical po- 
tentials of the reservoirs [|l^. Generalizing this method 
we calculate the spectral density of the two-level dot and 
the differential conductance for the current through the 
dot. As novel result we find a multiply peaked resonant 
structure both in the spectral functions and the differen- 
tial conductance. When the two levels in the dot are po- 
sitioned below the electrochemical potentials of the leads 
a triple-peak structure emerges in the nonlinear conduc- 
tance: a zero-bias maximum with one satellite peak on 
each side. When the levels lie above the electrochemi- 
cal potentials the peaks turn into notches. Depending 
on the level spacing our model interpolates between an 
5 = 1/2 and an 5 = 3/2 Kondo-type model (generahzed 
to include transitions between all 'spin states'). 

We denote the two energy levels of the dot by e+ and 
e_ (the choice of -I- and — arises from the eigenstates of 
the double dot system). Both levels are spin-degenerate, 
e^l = £il = £i for i S {-f , — }. The Hamiltonian is a gen- 
eralization of the Anderson model, Ti. — Tio.D + ^^o,^:cs + 

Ht- The term Ho,b = J2ia ^i<^4aC'ia + C^o"d("d - 1) 
describes the isolated quantum dot. The interaction of 
the electrons in the dot is accounted for by C/q, with 
?iD = J2i(7 '^ia'^ia being the electron number. The elec- 
trons in the reservoirs are taken to be non-interacting, 
Wo,ros = Eae{R,L}Efc<T£"fe'^4fea«afea- The tunneling 



Hamiltonian Ut = Eo.k^a{Tk^Aa"'o.ka + ^-c) describes 
tunneling between the dot and the reservoirs. 

The double-dot system coupled in series between two 
reservoirs as shown in Fig. Q with one spin-degenerate 
level in each dot is described by the Hamiltonian H = 
"^o.D + Wo, ICS + Ht + TIt- The Hamiltonian of the two 
capacitively coupled dots is Hq.d = J2da ^daC^i^c^^ + 
[71271171]. -I- C/o "'dT'^di- addition to the terms re- 
ferring to both dots it contains an additional interaction 
term Ui2n\n^. The index d € {1, r} denotes the left or 
right dot, respectively. The term Ht — So-l^^L'^i'a + 
h.c.) describes the tunneling between the dots, while the 
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tunneling Hamiltonian Wt is modified such that it only 
couples each dot to the lead adjacent to it. 

The double-dot model can be mapped onto the two- 
level single-dot model in the regime of large charging 
energy, i.e., when Uo and U12 are larger than all other 
energy scales of the system. Then it is sufhcient to 
consider only two adjacent charge states of the dou- 
ble dot. For appropriate values of a gate voltage cou- 
pled to the dots these states are the empty and the 
singly charged dot system. Then, we first solve the ex- 
act eigenstates of the double dot isolated from the leads, 
i.e., the eigenstates of -ffo,D + Hf. They are symmetric 
and antisymmetric combination of the single-dot states 
in the left and right-hand dots |-|-) = a\l) + f3\r) and 
|— ) = — a\r), with energies £+ and £_, respectively. 
The tunneling Ht couples the leads to these new eigen- 
states. Thus the effective tunneling matrix elements are 
Tl^+ = aT , rR,^+ = PT , Tl,- = (3T , Tr,,_ = -aT. 
Except for this modification of the matrix elements the 
double-dot system coincides with the two-level single dot. 

The mapping is not possible if we consider doubly- 
occupied states. The interaction between electrons in the 
same dot is stronger than between electrons in different 
dots. Hence two electrons in a double dot prefer to stay in 
different dots, and the low energy two-electron states are 
formed from these single-electron states only. In contrast 
in a single dot the two-electron states are composed of all 
single-electron states. By symmetry, in the case of three 
and four electron states the mapping is possible again. 

Physical quantities of interest can be expressed as 
quantum statistical averages. In particular, the d.c. cur- 
rent through the quantum dot takes the form 

I = Il = -Ir - - H.C.) . (1) 

The statistical average of an arbitrary operator 0{t) can, 
quite generally, be expressed as 



{0{t)) =tr\pQTKe^p 



-I j^dt'nT{t')i]0{t)i}. (2) 



The time integration is performed along the Keldysh con- 
tour, i.e., forward in time from some initial time U to t 
and then backward from t to U. The time-ordering oper- 
ator Tk arranges operators with respect to this Keldysh 
contour. The subscript I denotes the interaction picture 
with respect to Hq.d + 'Wo, res- 

The leads are assumed to be large and in thermal equi- 
librium. Their electrochemical potentials are kept fixed 
at /Xj. = — eVr (for r=L,R), and their electron distribu- 
tions are described with Fermi functions. We are inter- 
ested in the time evolution of the quantum states of the 
dot system, and therefore consider the reduced density 
matrix of the dot. To do this we first expand the expo- 
nential part in Eq. (^) into a series of the form 



Tk exp (^-i J dt' HT{t')?j = 



E 



(3) 
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Here the times of the integrals are ordered with respect to 
the Keldysh contour, such that ti > t2 > ■ ■ ■ > <m- The 
expectation value of this series is taken by tracing over 
the non-interacting reservoirs in the density matrices. We 
can perform the trace over the reservoir degrees of free- 
dom exactly since the unperturbed Hamiltonian Tio.ios is 
bilinear in the lead-electron operators. These facts allow 
us to apply Wick's theorem, and the forward and back- 
ward propagators become coupled due to the trace over 
the leads. In Refs. [^| Jl5|j2^ a systematic diagrammatic 
technique has been developed for presenting the terms 
of the expansion Eq. (||). There are well-defined rules 
for evaluating an arbitrary diagram and it is possible to 
identify and study processes of different orders. 

The current can be expressed in terms of the spectral 
densities Ai„{i^) of the electron states in the dot as |p3| 
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The factors F", = 2n J2k l^kral ^i^ - Eafea) describe the 
coupling between the dot level i and the lead a and they 
are assumed to be independent of energy Q . The spec- 
tral function describes the spectrum of possible excitation 
of the system. In general there exist separate Green's and 
spectral functions for the four states of the dot: one for 
each spin of each level. Here we restrict ourselves to the 
zero-field case such that the levels and the corresponding 
spectral functions are spin-degenerate. 

In general, we can not sum up all possible contribution 
from the expansion Eq. (^. Similar as in Ref. 
we proceed in a conserving approximation, taking into ac- 
count non-diagonal matrix elements of the total density 
matrix up to the difference of one electron-hole pair ex- 
citation in the leads. In this case we find a set of integral 
equations for the spectral densities Ai„{uj). While the 
coupling between the dot and the reservoirs is described 
within this approximation scheme, the double-dot states 
are solved exactly. This is reasonable when the coupling 
strength between the dots is at least comparable to the 
coupling between the dots and the leads. 

In the following, we choose a symmetrical bias /il = 
^/^R ~ eV/2 and consider equal-strength coupling of the 
dot to the leads Fl = Fr. The resulting conductance 
curves are symmetrical as a function of the voltage. The 
values for the energy levels that we mention, e.g., in the 
context of the illustrations, refer to the level positions 
before the renormalization due to tunneling. However, 
when we talk about levels being above or below the elec- 
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trochemical potentials we refer to the renormalized mea- 
surable level positions. Finally, we choose the temper- 
ature /cbT = 0.02r, where T Tl + Tr provides the 
energy scale, such that the thermal energies are always 
lower than the level spacing in the dots and much lower 
than the charging energy. 

In the spectral densities Ai(uj) there appear a com- 
plex, energy-dependent self-energy term which yields a 
broadening as well as renormalization of the dot levels. 
The energy-dependence of the self-energy term is due to 
the resonant-tunneling processes, and it results in a non- 
trivial functional form of the spectral functions. On top 
of this structure Kondo-type resonances may appear as a 
sharp peaked structure at the position(s) of the electro- 
chemical potential(s). These general effects of resonant 
tunneling have been found already in the investigations 
of a single one-level quantum dot |p^ , |l7| . 

We keep the lower level below the electrochemical po- 
tentials of the leads such that e+ ^ — F. Tuning the 
upper level e_ — e+ + Ae allows a crossover from the 
one-level limit, — > +00, to a fourfold-degenerate level, 
£+ = £_. For the single-dot model the level separation 
Ae reflects the actual physical level separation, whereas 
it is proportional to the coupling strength between the 
two dots for the double-dot model. In the limit of in- 
finite level separation the spectral density A^{lu) repro- 
duces results obtained before for a single level p^ , p7t . 
These are the level renormalization and broadening, as 
well as - for e+ ^ — F - Kondo-type resonances at the 
electrochemical potentials of the leads. For a low lying 
level this corresponds to the S = 1/2 Kondo problem, 
however generalized to a nonequilibrium situation if a 
transport voltage is applied. 

When the second level is included, the Kondo peak of 
the single-level model persists. However, the renormal- 
ization effect is stronger and, as long as £+,£- ^T, 
a further peak in the spectral functions appears (see 
the dashed curves in Fig. ||). This peak is shifted from 
the electrochemical potentials by Ae below (above) the 
Kondo peak in the spectral function for the lower (up- 
per) level. In the other extreme of zero level separation 
we have a fourfold degenerate case, and the model corre- 
sponds approximately to a spin S — 3/2 single-channel 
Kondo model (in the limit of low lying levels, further- 
more transitions between all 'spin states' are allowed in 
the present case) . In this case the resulting spectral func- 
tions are equal to each other and have a much more pro- 
nounced Kondo peak (see the solid curve in Fig. ^) than 
in the single level case. 

The large resonant peak, seen when = can be 
viewed as two Kondo peaks being on top of each other. 
When the degeneracy is lifted the peak splits into two, 
separated by Ae. One peak remains pinned at zero since 
the spin-degeneracy of the levels continues to lead to a 
Kondo peak. This result differs from the situation where 
the spin degeneracy of a single level is lifted by a magnetic 



field |T|l. When the upper level is lifted far above the 
lower one the properties of the system approach those of 
a single-level dot thus exhibiting just one Kondo peak. 

The spectral functions are normalized. Since we have 
neglected multiple occupancy, the normalization is 

/oo 
duj Aia{i^) = PO +Pia , (5) 
-00 

where po and pia- are the probabilities of an empty dot 
and that with single electron state ia occupied. This 
expression allows us to relate differences in the magni- 
tudes of the spectral functions to differences in occupa- 
tion probabilities, e.g. when the spectral functions corre- 
sponding to different levels are compared to each other. 

The bias voltage has the effect of splitting all the res- 
onant peaks in the spectral functions by eV, see Fig. ||. 
The peaked structure of the spectral functions is reflected 
in a novel structure in the nonlinear conductance, see 
Fig. IJ. Instead of a single zero-bias peak, known from 
the single level dot, there appear now three sharp res- 
onant peaks between the broad peaks corresponding to 
the renormalized energy levels. The zero-bias anomaly, 
characteristic for the single level case remains, but two 
additional satellite peaks appear at eV = ±Ae. When 
the level separation is increased the new peaks decrease 
in size while the peak at zero bias remains unchanged. 
When the levels lie above the electrochemical potentials, 
the nonlinear conductance exhibits a triplet of minima 
with one minimum at the electrochemical potentials and 
two satellites on each side. When the levels are on op- 
posite sides of the electrochemical potentials the spectral 
functions as well as the conductance show very varied 
structures, which are not easily classified. 

The systems considered here have not yet been inves- 
tigated experimentally. Ralph and Buhrman have 
measured tunneling ~ enhanced by the Kondo effect - 
through a localized state of an impurity equivalent to a 
single one-level dot. Some unexplained features in their 
results may be due to the existence of the other levels 
above the electrochemical potentials. The model of a 
single- or two-level dot with only one electron in it could 
possibly be realized with the setup used by Tarucha et 
al. . They argue that they can resolve situations with 
one or few one spin-degenerate levels and one or few elec- 
trons in the dot. In particular, they report that a third 
electron inserted into the dot finds two spin-degenerate 
levels whose separation can be continuously varied by a 
magnetic field. In this system, if the gate voltage is fixed 
near the first or third conductance peak the single- or 
double-level dot is realized, respectively. The difference 
should be observable in the nonlinear conductance pro- 
vided that the zero-bias anomalies can be resolved. The 
relevant parameters reported in Ref. pO[ | - the coupling 
strength (tunneling rate) F, voltage, and temperature - 
are already close to those used in our calculations. 



3 



In conclusion, we have described electronic transport 
through a two-level quantum dot and a single-level dou- 
ble dot coupled between reservoirs. A real-time diagram- 
matic technique provides a very general tool for the de- 
scription of various systems, both in the perturbative and 
nonperturbative regime. We have found a novel triple- 
peak structure in the nonlinear conductance of these sys- 
tems and suggest experiments where it is observable. 
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FIG. 2. Spectral function A_|_(o.)(a;) with £+ = —2. OP and 
£- = e+-|-Ae. The large Kondo peak of the four-fold degener- 
ate level splits when the degeneracy is lifted. The separation 
equals the level spacing, Ae e {0, 0.2r, 0.4r}. 

FIG. 3. Spectral functions A+{uj) when a bias voltage V is 

applied. The voltage splits all the resonant peaks by energy 
eV and the peaks cross each other when eV = Ae. 

FIG. 4. The figure shows apart from the broad main peak 
in the non-linear conductance, which is sometimes referred to 
as the resonant peak, sharper resonant peaks close to zero. 
The inset shows the same enlarged. The solid curve is for the 
four-fold degenerate level. The dashed curves are for split lev- 
els, Ae = 0.2r (long dashed) and Ae = 0.4r (short dashed). 
The conductance g{V) is expressed in units of e^/h. 
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